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Abstract 

j^-^ ' We consider lattice equations on 1? which are autonomous, afBne linear and possess the sym- 

, metries of the square. Some basic properties of equations of this type are derived, as well as a 

sufficient linearization condition and a conservation law. A systematic analysis of the Lie point 
and the generalized three- and five-point symmetries is presented. It leads to the generic form of 
, the symmetry generators of all the equations in this class, which satisfy a certain non-degeneracy 

condition. Finally, symmetry reductions of certain lattice equations to discrete analogues of the 
Painleve equations are considered. 

1 Introduction 
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J> ' The importance of symmetry based techniques appUed to differential equations, especially to nonlinear 

, ones, is well known. It can be argued that symmetry methods are the most effective ones for obtaining 

explicit solutions of complicated (systems of) nonlinear partial differential equations. In fact, the so- 
called group invariant solutions of such equations form a well known example of the results arrived at 
by these methods. It is also well known that, using its symmetries one can construct new interesting 
solutions of a given equation from much simpler ones. 
' However, the notion of a symmetry group of transformations acting on the solution space does not 

, have to be limited only to differential equations. It can be equally well applied to other types, such 

as the algebraic and difference ones. The latter arise in many diverse branches of mathematics and 
physics, such as discrete geometry, integrable systems, special functions and orthogonal polynomials, 
the study of exactly solvable models in statistical mechanics, crystal lattice theory and many others. 
The wide range of their applications shows that difference equations are of equal importance with their 
differential counterparts. As a result, symmetry methods have started being applied to the analysis of 
difference equations, as well. 

In fact, difference equations have already been studied by symmetry methods from various points 
of view, see e.g. O [7l [3 [TU [TBI [El [13 [26] and references therein. Symmetries of integrable partial 
difference equations first appeared as compatible constraints in the work of Nijhoff and Papageorgiou 
[20] . The motivation was a specific reduction of the discrete modified Korteweg-de Vries (KdV) equation 
to a discrete analog of Painleve II, in the same spirit as the Painleve II ordinary differential equation 
arises as a similarity reduction of the celebrated partial differential equation of KdV. Further examples of 
such compatible constraints for integrable partial difference equations were given in [TUl [IHl [SB [ISl [30] • 
In this paper, we present a systematic study on the symmetries and reductions of autonomous partial 
difference equations, which are not necessarily integrable. Specifically, we consider a quite large class of 
lattice equations defined on an elementary quadrilateral, which contains the integrable ones classified 
recently by Adler, Bobenko and Suris in [1]. The members of this class are characterized by (i) affine 
linearity and (ii) D4-symmetry i.e. the symmetries of the square. 
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We first prove that each equation in this class admits two, at least, three-point generalized symme- 
tries. They are determined by a pair of polynomials arising from the equations defining relation. For 
the generic case, we give the form of the generators of three- and five- point generalized symmetries, 
and a greatly simplified form of the corresponding determining equations. These results also extended 
to symmetry transformations which act on the lattice parameters appearing in the equation, as well. 

The previous symmetry analysis is then applied to the equations obtained in classification [1]. The 
result of this part of our investigation is an exhaustive list of the corresponding Lie point, three- and five- 
point generalized symmetry generators. As a final application, we consider specific symmetry reductions 
of the discrete potential KdV to ordinary difference equations, which represent discrete analogues of 
those of Painleve. 

The paper is organized as follows. Section [2] contains the necessary preliminaries on symmetries of 
difference equations and the notation that we use in the following sections. In Section [3] we introduce 
the family of lattice equations under consideration and its characteristic properties. The main results 
of the symmetry analysis are presented in the following four sections, where the general form of the 
corresponding symmetry generators is given explicitly. The Lie point symmetries are studied in Section 
m The three- and five-point generalized symmetries are presented in Sections [5] and [H respectively. 
In Section [7] we extend the previous considerations to symmetry transformations acting on the lattice 
parameters, as well. In Section[8]we present the symmetries of the equations of the classification [T], and 
Section [S] deals with symmetry reductions. We conclude with Section [TUl where an overall evaluation of 
the results obtained in the main body of the paper is presented, along with various perspectives on the 
subject. In the Appendix, a detailed proof of the Proposition of Section [5] is given. 

2 Preliminaries on symmetries of difference equations 

A partial difference equation is a functional relation among the values of a function u : Z x Z ^ C (or 
CP) at different points of the lattice, which in general involves the independent variables n, m and the 
lattice spacings a, /3, as well, i.e. a relation of the form 

£" (n, m, M(ri, m), u(n -f 1, m), m -f 1), . . . ; a, /3) = 0. (1) 

The analysis of partial difference equations is facilitated by the use of two translation operators on 
functions on 1? , defined by 

{^ll'^uj (n, m) = u{n + k,m) , (^S^^^uj (n, m) — u{n, m + k) , where G Z . 

It is also found useful to introduce the notation 

u{n, m) — M(o,o) : 'fJ'ik.i) — u{n + k,m + I) , where k, I G , (2) 

for the values of the function u and this will be adopted from now on. 

Let G be a one-parameter group of transformations acting on CP, the domain of the dependent 
variable U{o,o) of a lattice equation, i.e. 

G : U(o,o) ^'(", W(o,o);e) , e £ C. 

We denote by J'^'^) the forward lattice jet space of order fc e IN with coordinates j)), where i,j G IN 
and i + j < k . Similarly, one can define the backward lattice jet space of order fc, denoted by J'^"'^-', 
with coordinates i,j e and i + j < k, and in general the fc-order lattice jet space 

with coordinates (u(±i,±j)), i,j G IN and i + j < k. The prolongation of the group action of G on J'^'^-' 
is defined by 

G^'^') : ^ ($(n + i,m + j,'U(,_j);e)) . (3) 
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The infinitesimal generator of the group action of G on the domain of the dependent variable is given 
by the vector field 

X = i?(n,m,U(o,o))9„(„ „, , 
where the symmetry characteristic R{n,m,U(^Q^o)) is defined by 

The group action is reconstructed by exponentiating that of the vector field x 

$(n,m,U(o^o);e) = exp(ex) M(o,o) • 

The infinitesimal generator of the action of G'^'^' on J*^*^' is the associated k^^ order forward prolonged 
vector field 

fe fe- 



X 

1=0 j=0 



J2 E {^n ° S^^R) {n, m, «(o,0)) 



The transformation group G is a Lie-point symmetry of the lattice equation |T]) if it transforms any 
solution of ^ to another solution of the same equation. Equivalently, G is a symmetry of equation ([T]), 
if the latter is not affected by the transformation ([3|). The infinitesimal criterion for a connected group 
of transformation G to be a symmetry of equation ([1]) is 



(n,TO, M(o,o),'«(i,o),M(o,i), • ■ ■)) =0- (4) 



This should hold for all solutions of equation ([!]) and, thus, the latter and its consequences should be 
taken into account. Equation ^ delivers the most general infinitesimal Lie point symmetry of equation 
P^. The resulting set of infinitesimal generators forms a Lie algebra g from which the corresponding 
symmetry group G can be constructed by exponentiation. 

A lattice invariant under the action of G is a function / : — > C which satisfies /(g'-'^-' • 
("(±i,±i))) = I{u(±i,±j)) for all 51 G G and all (M(±i,±j)) G jC^'.-'^'). For connected groups of transforma- 
tions, a necessary and sufficient condition for a function / to be invariant under the action of G is the 
annihilation of / by all prolonged infinitesimal generators, i.e. 

= 0, (5) 

for all X G g. 

By relaxing the geometric assumption that the symmetry characteristic R depends on n, m and 
U(o,o)j only, and allowing i? to be a function defined on Z2 X for some finite but unspecified 

/c G A: > 1, we arrive naturally at the notion of the generalized Lie symmetry. Symmetry generators 
of this type cannot be associated to transformation groups acting geometrically on the domain of the 
dependent variable. Lowest order (k = 1) generalized symmetries are given by the following vector field 



V = -R(n,TO,U(o,o),'U(i,o),'"(o,i)''"(-i,o)iW(o,-i))c^«(, 



0,0) 



3 A class of two-dimensional lattice equations 

In this section, we present a class of two-dimensional lattice equations, which involve the values of a 
function u at the vertices of an elementary quadrilateral as shown in Figure [1] Specifically, we consider 
the two-dimensional lattice equations of the form 

<3("(0,0)>W(l,0):'"(0,l)''"(l,l);'^'/3) = 0' (6) 

where the function Q 
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"(CO) a M(i,o) 
Figure 1: An elementary quadrilateral 



• does not depend explicitly on the discrete variables n, m, 

• depends explicitly on the values of the unknown function u at the vertices of an elementary 
quadrilateral, i.e. c^«(i jj Q(w(o,o)) '"(i,o)> ^^(0,1)1 ""(i.i); ctj /3) 7^ 0, where i, j = 0, 1, and may depend 
on the parameters a, (3 of the lattice, 

• is linear in each argument (afhne linear): 9^^. (5(m(o 0)7 0) ! W(o 1), U(i i); a, /3) = 0, where i, j = 
0,1, 

• and possesses the symmetries of the square (D4-symmetry) : 

(3(^^(0, o),W(i,o),M(o,i)'"(i,i);"'/^) = eQ(M(o,o), "(0,1)' ^(1,0)7 "(1,1);/^'") 

= Cr(3(M(i,o), "(0,0), "(1,1)' "(0,1)! P) ■> 

where e = ±1 and a = ±1. 

The symmetry analysis of the above class of equations is significantly simplified by the use of certain 
polynomials arising from the function Q. In the rest of this section, we define these polynomials and 
derive some of their properties, in order to make the symmetry analysis of the following sections more 
concise. 

To begin with, we note that the linearity of the function Q implies that the functions 



^("(0,0), "(1,0); Q^, /?) = Q Q,«(o,i)«(i,i) - '9,«(o,i)<3,M(i,i) ' (7a) 

^i("(o,o): "(0,1); a, /?) = Q Q,«(i,o)"(i,i) ~ Q.uci.ojQ.fci.i) ' ('''b) 

^2("(o,i), "(1,1); a, /?) = Q Q,«(o,o)«(i,o) ~ '5,«(o,o)<3,«(i.o) 5 C^c) 

'i3("(i,o), "(1,1); a, /?) = Q Q,«(o,o)«(o,i) ~ '5,«(o,o)<3,«(o.i) 5 ('''d) 

are biquadratic polynomials in their two first indicated arguments, and the same holds for the functions 

^("(0,0)) "(1,1); Q^, /^) ~ <3Q,tl(i,0)M(0,l) ~ ^."(l.OjQ.MCO,!) 5 (8a) 

G'i("(i,o),"(o,i);a,/3) = <3Q,«(o,o)"(i,i) ~ ^."(o.ojQ.uci.i) • (8b) 



In fact, the linearity of the function Q and the above definitions lead immediately to the properties 
expressed by the following two lemmas. 

Lemma 1 Let the Junction Q be ajfine linear. The polynomials defined by are constants if and 

only if the function Q is linear, i.e. 

Q = fi{a, /3)"(o,o) + /2(a, /3)u(i^o) + fsia, /3)"(o,i) + /4(a, /3)"(i,i) + fsi^, P) ■ (9) 

Proof : If the function Q is of the form Q then definitions ([7]), ^ imply that these polynomials are 
constants. Conversely, assuming that these polynomials are constants, we solve (O, (|S]) for the second 
order derivatives of Q and take the compatibility conditions among the resulting equations. This leads 
to eight first-order partial differential equations for the function Q. This overdetermined system of 
partial differential equations imply that Q is necessarily of the form □ 
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Lemma 2 // the function Q is affine linear, then the relations 

/i(w(o,o): "(1,0); a, P) /i2("(i,o), ""(1,1); a, P) = ^i(w(o,o), "(0,1); P) ^3(^(1,0), ""(i.i); P) 

= G(w(o,o),U(i,i);a,/3)Gi(u(i_o),W(o,i);a,/3) (10) 

hold, in view of the egttaiion (5(u(o.o): ^(1.0)1 W(o. 1)1 '^'(1,1); c^? /3) =0. 

Proof : It follows from the definitions of the functions involved, the affine linearity of the function Q 
and by taking into account the equation Q = 0. □ 
The D4-symmetry of Q implies that the polynomials defined by ([Tj) are symmetric in their first two 
arguments and related as follows 

'i("(o,o),'"(i,o);"7/3) = /ii('"(o,o),'"(i,o);/3, a) = ^2(u(o,o), W(i,o); /?) = ^3("(o,o), "(1,0); (H) 

Moreover, the functions G and Gi, defined by ([8]), are symmetric in their first pair of arguments and in 
(q!,/3), and they have the same form, i.e. 

Gi(w(i,o),U(o,i);Q^'/5) = G'(w(i,o),M(o,i);a,/3) ■ (12) 
In this case, Equation (fTO]) simplifies to 



/l("(0,0),U(l,0)) ^("(0,1), "(14)) = /l("(0,0),'«(0,l)) ^("(1,0), "(1,1)) 

= 0(^(0,0), u(i4))G'(-U(i^o), "(04))- (13) 

Remark 1 In the following, we omit the dependence of the polynomials h, G on the lattice parameters. 
It should be noted that when the polynomial h is evaluated at two neighboring points in the horizontal 
direction of the lattice, then the parameter dependence is (ck,/3) and the order is reversed when h involves 
two points in the vertical direction. 

Remark 2 // the function Q is affine linear and possesses the Hji- symmetry, then, following the proof 
of Lemma 2, one can prove the validity of the following relations 

Q\ ^ M^(i.o)^j;(M))gKi,o).mi)) ^ (14^) 
^ Mmo).mi))gKi.o).mi)) ^ (^4^) 



?,"(i,o) ^(W(0,1),U(1,1)) 



G(M(i,o), U(0,1)) 
<3,"(0.i) ^1(^(0,0), U(1,0)) 



0, (14c) 
0. (14d) 



Q,«(0,0) G'(M(i,o);"(04)) 

The latter are quite useful in the symmetry analysis of the equations under consideration. 

Relations hold, in general, in view of the equation (5 = 0. However, in certain cases, these 
relations hold identically, i.e. without taking into account the equation Q ~ Q. In such cases, the 
corresponding equations can be linearized using an appropriate transformation. 

Proposition 1 Let the function Q be affine linear and possess the Y>i-symmetry. If the relation 

^("(o,o),M(i,o))^("(o,i)''"(i,i)) ~ '^('"(o,o),w(i,i))G(m(i,o),U(o,i)) = (15) 

holds identically, i.e. without taking into account the equation Q — 0, then the polynomials h, G are 
factorized as 

^("(o,o):'«(i,o);a,/3) = p(M(o,o);a,/?)p(w(i,o);a,/3) , 
G(u(i,o),M(o,i);a,/3) = ±p(u(i,o);a,/3)p('"(o,i);a,/3) , 
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and the equation Q — is transformed to a linear equation under the transformation 

f 1 

u — > u — T(u) :— / — TK^u. 

J p{u;a,l3) 

Proof : First, we write equation (fT5|) in the form 

fe(M(0,0)7'»(l,0)) _ G'(m(o,o),M(i4)) 

G(w(i,o),U(oa)) /i(M(o,i)''"(i,i)) ' 

Differentiating this relation w.r.t. U(i.o) (equivalently w.r.t. M(i i)), we find that the polynomials h, G 
must have the form 

h{'^(ofi),U(ifi)]a,l3) = p(u(o,o);a,/?)p(w(i,o);a,/3) , 
G(u(i^o),'"(o,i);a,/5) = ±p(u(i^o);a,/^)p(w(o,i);a,/3) , 

where p is a quadratic polynomial in its first argument and symmetric in (a, /?). 
Now let ^ 

u = r(u) := / — ^du, 

and 

— 1 /2 

-F(M(o,o),"(i,o),"(oa),w(i,i)) = Q • (p(^i(o.o);a,/^)p(w(i,o);a,/3)p(u(o4);"'/5)p("(i,i);"'/^)) 

Using the previous result, it immediately follows that, under the transformation {u,Q) — *■ (u,F), the 
partial differential equations (O, ^ simplify to the following system 



FF - 


~ ^,"(0,1) -^."(1,1) 


= 1, 


FF - 


~ ^i"(l,0) -^."(1,1) 


= 1, 


^ -^."(0.0)"(1,0) 


~ ^,"(0.0) -^:«(1,0) 


= 1, 


^ -^."(0,0)"(0.1) 


~ ^,"(0,0) -^.«(o,i) 


= 1, 


^ ^ :"(1,0)"(0,1) 


~ ^,«(1,0) -^."(o.i) 


= ±1, 


^ >"(0.0)"(1,1) 


~ -^>«(0,0) ^,"(1,1) 


= ±1 . 



From this it follows that the equation -F(u(o.o); '^(i.o)? 'I'co,!)! — must be linear. C 

Example 1 Consider the equation 

("(i,o)"(04) + '"(i,o)W(i,i) + ""(0,1)^^(1,1)) + "(i,o)W(o,i)W(i,i) + "(0,0) + "(1,0) + U(o,i) + w(i,i) = 0, 
listed in [T3]. One easily finds that, in this case, 

^('^(o,o),'«(i,o)) = -("(0,0) - l)Ki,o) - 1) ' G(m(i,o),U(o,i)) -("(1,0) - l)Ko,i) - 1) • 



Preforming the transformation 

u 

or, equivalently, 

the above equation linearizes to 



du 1 , f u — 1 



f 2 



"(0,0) + "(1,0) + "(0,1) + "(1,1) 
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On the other hand, if the relations (fT3l do not hold identically, then the equation Q = can be 
written as a conservation law. Specifically, 

Proposition 2 Let the function Q be affine linear and possess the D 4- symmetry. If the relation 

^("(0,0), ■"(1,0)) ^("(0,1), "(1,1)) = ^("(o,o),W(o,i)) ft.(u(i^o),"(i,i)) (16) 

does not hold identically, then the equation Q — can be written in the form of a non-trivial conservation 
law 

(5,„ - irf) Fi(n,TO,'U(o,o),U(i,o)) = (^n - ^2(71, m, U(o^o), "(0,1)) , (17) 

where 

Fi(n,m,U(o^o),W(i,o)) = (-1)"^" In /i(m(o,o) , "(1,0)) , (18a) 
F2(n,m,U(o^o),U(o,i)) = (-1)"^" In /i(u(o,o) , "(ci)) • (18b) 



Proof : It follows by combining (fT7| with ([T8|) and taking into account relation (fT6|) . □ 

4 Lie point symmetries 

In this section we study the Lie point symmetries of the two-dimensional lattice equations under con- 
sideration. Let 

X = (/)(n,TO,M(o,o))9„(„„, 
be the generator of a point symmetry transformation of the equation 

<9("(o,o),M(i,o),'«(o,i),M(i,i);a,^) = 0. 
The infinitesimal symmetry criterion 

x^^) (Q(u(o,o),'"(i,o),U(o,i),w(i,i);a,/3))|Q^o " ^' 

where 

1 1 

x(2) ^ ^^(P{n + i,m + j, U(, ; a, (3) du^,.^ , 

implies that the determining equation 

'3."(o,o)<^("'"^'"(o,o);a7/3) + Q,«(i,o) '/'('■* + 1'"^' "(1,0);": /3) 

+ <3,"(o.i)<^("-''" + 1'"(o,i);q^;/5) + Q,«(i,i)<^(" + 1>"^ + i>"(i,i);";/5) = (i9) 

should hold on every solution of the equation Q — 0. 

Since the function Q is linear in U(^i i-^, the equation Q = can be uniquely solved for it^j^ x) in terms 
of W(o.o), W(i.o) and U(o.i)- Using the relations 



<9,"(i.o) _ ^'("(0,0), "(0,1)) Q,u^o,i) _ ^("(0,0)1 "(1,0)) '9,"(i.i) _ '^."(1,1) 



(20) 



Q,^{0.0) G'(m(i^o),"(0,1))' Q,«(o,o) G(u(i^o),'"(0,l)) ' Q,«(o,o) G(li(i^o),"(0,l)) ' 

we eliminate it (1,1) from equation (|19p and arrive at 

J.( »\ , '^(■"(0,0), "(0,1)) , -, m 

0(n,m,U(o o);a,p) + 777 -(/)(n 1, m, M(i q); «, P) 

G(u(i^o),"(o,i)) 

, ^("(0,0), "(1,0)) , , J,^ I 1 11 m ^01^ 

+ 777 r<?'("-,?7T- + l,"(o,i);a,p) = 7w r0(n + l,m+l,U(ii);a,/^). (21) 

G("(i,o),"(o,i)) G(m(i,o),"(o,i) 
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Since h is, in general, a quadratic polynomial of U(o,o)j differentiating the determining equation (|2ip 
three times w.r.t. U(o,o) smd once w.r.t. it(i_o), we arrive at 



where Dj denotes the total derivative operator, i.e. 

Dj = dj + , where j = ^(0,0), "(1,0), "(o,!) • 

Writing equation explicitly, one arrives at 

^«(i,i) (^("(oa),"(ia))^5u(i,i) + 1,"^ + =0- (23) 
The last equation splits into the following system 

^«<o,i,(0,W(i,i))<,i,<^ + 2^«,o,i,«a,i,(0,W(i,i))5«a,i)'^ = ' (^4) 

^,«(O.l)«(O,l)(0''"(l,l))^«(l,l)<?^ + 2/l,U(o,i)U(o.l)«(l,l)(0''«(l,l))^U(l.l)'^ = 

where we have omitted the arguments of the function (f>(n + 1, m + 1, a, (3). If the matrix 

/i(0,U(i,i)) i)(0,'U(i,i)) 
^ = I ^«(o,i)(0>w(ia)) ^,«(o,i)«(i,i)(0,"(i,i)) I (25) 

^:«(0,1)"(0,1) (f^' ""(1,1)) ''•■"(0,1)1'(0,1)'"(1,1) ('^' ""(1,1)) 



has rank equal to 2, then the system (p4)l has the unique solution 

+ 1, TO + 1, U(i 1); a, /?) = 9^(/)(n + 1, m + 1, a, /3) = 0, 

leading to 

m, U(o,o);a,/3) = v42(", "i; a, /3)m(o q-, + (n, to; a, /5)u(o,o) + Ao(n, to; a, /5) . 
In the degenerate case, where rank;B = 1, the function /i(u(o,i), U(i,i)) separates variables, i.e. 

/i(w(o,i)''"(i,i)) = 'io("(o,i))'*o(m(i,i)) , (26) 

where the function is, in general, a quadratic polynomial of its argument. In this case, the system 
reduces to one single equation, namely 

(^/io(u(i4))^9,^j^_^j0(n + 1,TO,+ l,U(i^i);a,/3)^ = 0, 

which integrated once yields 

^3 ,/ , T ,1 m ^3(71+ l,m+ l;a,/3) 

< • ' /io(m(i,i)) 

Once the function Q is given, the last equation can be easily solved leading to the general form of the 
corresponding symmetry characteristic. 

The form of the functions Ai(n,m;a, P) is obtained in the following way. We first substitute the 
resulting characteristic of the symmetry generator into equation (|2ip and then use the equation to 
eliminate U(^i iy Setting equal to zero the coefficients of the different monomials of the remaining 
variables, we finally arrive at a linear overdetermined system of difference equations for the unknown 
functions Ai(n,m;a,P). The solution of this system delivers the Lie point symmetries of the lattice 
equation Q = 0. 
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Remark 3 Actually, the above symmetry analysis can be immediately extended to the case where the 
function Q is affine linear but not D4 symmetric. The procedure followed above leads to, essentially, the 
same equations. They are the ones obtained from equations \21\] and V23\] by making the replacements 

Ku{id.,o),u(^i,Q)) ft,(u(o,o),M(i,o);a,/3) , /«(w(o,o) , "(0,1)) /ii("(o,o), "(0,1); /3) 1 
G'(u(i,o),'«(o,i)) ^ Gi(u(i^o),'«(o,i);a,/3) , 



and 

respectively. 



'^('"(0,1)' "(1,1)) ^2(m(o,i)''"(i,i);"'/5) 



5 Three-point generalized symmetries 

In this section, we consider more general symmetries than the point symmetries presented in the 
preceding section. More precisely, we search for generalized symmetries with characteristics of the 
form i?(n, m, U(o,o): ■'^(1, 0)1 "(-1,0)! Q^j Once the symmetries of this type are found, one can con- 
struct similar symmetry characteristics in the other direction of the lattice, simply by interchanging 
mutually the lattice variables and parameters. This follows from the D4-symmetry of the function 
Q('"(o,o)j ""(1,0): "(0,1): 1); a, (3). Thus, it suffices to restrict our considerations to symmetry character- 
istics of the form i?(n, m, M(o,o), M(i,o)5 W(_i^o); ct, see Figured For brevity, we refer to this type of 
symmetries as three-point generalized symmetries. 

^(-1,0) ^ (0,0) ^(1,0) 
Figure 2: Three consecutive horizontal lattice points 

The three-point generalized symmetry analysis for the equations under consideration is summarized 
in the following 

Proposition 3 Every two-dimensional lattice equation Q{u[q (1.0)1 P) — Oi where the 
function Q is affine linear and possesses the 1^4^- symmetry, admits a three-point generalized symmetry 
with generator 



"(1,0) - "(-1,0) 

Moreover, in the generic case, where the matrix 



^/i,«(i,o, (w(o,o) , W(i.o); a, 9„,, 



0) 



h{x,y) G{x,z) G{x,w) 
G = I h^x{x,y) G^x{x,z) G^x{x,w) 

h.xx{x,y) G,xx{x,z) G^xx{x,w) , 

has rank 3, every three-point generalized symmetry generator necessarily has the form 



(28) 



V„ = a(n; a, /?) v„ -f ^^(n, m, U(o^o) ; /3)5«(o,o) ' 

where the functions a{n] a, (3), (j){n,m,U(^Q Qy,a, (3) satisfy equation \29fl . below. 

The Hji-symmetry of the function Q implies that the vector fields resulting from the mutual replace- 
ments 

n < — > m , a < — > (3 , U(i^o) ' — > "(o,i) ■, 
in v„, V„, are the generators of the three-point symmetries in the vertical direction. 
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Since the proof of the Proposition is involved, we give here only a sketch of the proof and leave the 
details for the Appendix. 

1. We differentiate the determining equation w.r.t. to W(i _i) and M(i,o) ^nd find that an admissible 
symmetry characteristic has the form 

'^("(0,0), "(1,0);",/^) 1, I . /^^^ , <?^('^'™''"(o,o);W(i,o);a,/3) 



A(ri, m; a, /3) (u(o,o),W(i,o);a,/5) + 

V "(1,0) ~ "(-1,0) ^ ' / 



where A and (f) are arbitrary functions of their arguments. 

2. Next, we prove that the function A is independent of the discrete variable m and, consequently, 
the determining equation is independent of g) ^.nd i). 

3. Then, we show that the determining equation, apart from the function 0, is also independent of 

"(2,0)- 

4. Finally, we eliminate the value W(i,i) from the determining equation, using the equation Q = 0. 
At this point, wc conclude that every equation, which is afhne linear and possesses the symmetries 
of the square, admits a three-point generalized symmetry with generator 

//i(u(o,o),M(i,o);a,/?) 1, . m\ a 

= — :5/i,«(i,o)("(o,o),M(i,o);a,/3) a„ . 

\ "(1,0) - "(-1,0) ^ / 

5. In the generic case, where the matrix Q has rank 3, we prove that the characteristic of the symmetry 
generator R{n, m, M(o.o) , "(i.o) ■, "(-i.o) ]ol,(3) necessarily has the form 

t /^(■«(o,o),W(i,o);a,/3) 1, , ri^\ , (l>{n,m,U(^aoy,a,j3) 
R = a[n;a,(J) -/i,„ (u(o.o) , U(i,o); p) H 7; . 

\ "(1,0) — "(-1,0) ^ / ^ 

where the functions a{n;a,(3), 4>{n,m,U(^ofiy, a, P) satisfy the simplified determining equation 

w a\ / , T ht ^2^^ f ^("(1,0)7^(0,1)) 

\ "l"(0,0)i "(i,o)J 

+G(u(i^o),"(o,i))'^("-,'^, "(o,o);a,/3) + /i(u(o,o), "(o,i))</'("- + ^,m,U(i^Qy,a, /3) (29) 
+'i("(o,o),"(i,o))0(","7, + l,U(o,i);a,^) = Q'^^^^_^^(p{n + l,m + l,U(^i^iy, a, P) . 

6 Five-point generalized symmetries 

In this section we consider symmetries with characteristic depending on the values of u assigned on five 
vertices, which form a cross configuration on the lattice as shown in Figure [31 We refer to this type 
of symmetries as five-point generalized symmetries. Such symmetry generators exist since any linear 
combination of the form civ„ -|-C2Vm, where v„, v™ are given in the previous section, is also a symmetry 
generator. 
Let 

1 1 

,^^(2) = ^ ^ R{n + i,m + j, u^^ij) , j) , "(» j+i) , "(j-i,j) , "(j,j-i) ;«,/?) 

be the prolonged generator of a generalized symmetry of the equation Q — 0. Acting on the latter with 
w^'^-' we obtain the determining equation 

'9,"(o,o)^("' "(0,0), "(1,0): "(0,1), "(-1,0), "(0,-1); "7 P) 

+ '9,«(i,o)^(" + 1, "(1,0), "(2,0), "(1,1), "(0,0), "(1,-1); a, P) 
+ Q,«(o,i)^(»^, ™ + 1, "(0,1), "(1,1), "(0,2), "(-1,1), "(0,0); a, (3) 
+ <3,"(i,i)^(^i + 1,"^+ l,"(i,i),"(2,i),"(i,2),"(o,i),"(i,o);a,/3) = 0- (30) 
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■"(04) ^ 



"(0,0) 



•"(1,0) 



• • • 

"(0,-1) 

Figure 3: Cross configuration 

Since U(2,i) depends on U(2.o) through the equation Q(m(i,o)j "(2,o)j "(i,i)i "(2,1); c^j /?) = 0, and U(i 2) 
depends on M(o.2) through the equation Q(m(o,i)7 U(i 1), U(o,2)j "(1,2)) Q^i /?) = 0, differentiating equation 
(|30p once w.r.t. U(2,o) ^-nd then w.r.t. W(o,2)7 we arrive at 

dui2 1) C^"(l 2) 

-R,«(2i)«(i2)('^ + l'"i+l>"(i,i);"(2,i),"(i,2),'"(o,i),"(i.o);a,/3)^ — ^ ^ — ^ = 0. (31a) 

OU(2,0) C'"(0,2) 

On the other hand, 0) depends on i) through the equation Q('U(_i o)j "(0.0)5 "(-1,1) i "(0,1)! 0^1/?) = 
0, and U(o,-i) depends on U(^i _i) through the equation Q(u(-o.„i), M(o,o)5 "(i,o)i Q^i /5) — 0- Thus, 

the differentiation of equation ([50)1 once w.r.t. 1) and then w.r.t. U(i „i) yields 



duii 0) 9u(o -1) 

^,"(-1.0)"(0,^l) i'^^ "(0,0) : "(1,0): "(0,1): "(-1,0) 7 "(0,-1) ; P) g^^ ^ q^^^ = 0- (31b) 



Moreover, differentiating equation ([50]) once w.r.t. W(2,o) (respectively U(o.2)) and then w.r.t. U(i _i) 
(respectively j^)), we obtain 

^,«(2.o)«(i.-i)("- + l>"i'"(i,o),"(2,o),"(i,i),"(o,o),"(i,-i);a,/3) = 0, (31c) 
-^,"(o,2)«(-i.i)("'"^ + l'"(o,i)'"(i,i)'"(o,2),"(-i,i),"(o,o);a,/3) = 0- (31d) 

Equations ()31|) imply that 



i?(n,m,M(o,o),U(i,o),W(o,i),"(-i,o),"(o,-i);a,/3) = ■Ri(", U(o^o), "(i,o), "(-i,o); /3) + 

i?2(n, m, U(o,o): "(0,1) > "(o,-i) ;«,/?)■ 

Substituting the above relation into the determining equation pop . and following the steps of the proof 
in Section [5l we find that the function Ri has the form 

u( R\ A( q^ fe("(o,o),"(i,o);«,/?) , ^ t r,^ 
Ri (n, m, U(o,o) , "(i.o) , "(-i,o) ; P) = A(n, m;a,fi) h (pi (n, m, U(o,o) , W(i,o) ; «, p) , 

"(1,0) - "(-1,0) 

where A, are arbitrary functions of their arguments. In a similar manner, one finds that 

h{u/Q Q),urQ i)] [3, a) 

i?2(n,m,U(o,o),M(o,i),"(o,-i);a,/5) = B{n,m;a,P) ■ ■ h <?i2(n, m, -(/(o.o), "(o.i); /?) , 

"(0,1) — "(0,-1) 

where B, (j)2 are arbitrary functions of their arguments. 
Setting 



= -^A(n,m;a,/3)/i,„(^ Qj(u(o,o),M(i,o);a,/3) + ^ '/'i(", "t,, U(o,o), U(i,o); /?) 
1 , , ^ , 1 



02 = - - B(n,m;a,/3) /i,tj(o j)(u(o,o),M(o,i);/3,a) + - V'2(^^, m, U(o,o), U(o,i); /3) : 
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we substitute the above relations into the determining equation and differentiate the result w.r.t. M(2,o)- 
This leads to the constraint 

A{n + 1, m; a, (3) — A{n + 1, m + 1; a, /?) = , 

which implies that 

Ain, m; a, P) = a{n; a, /3) . 
Similarly, differentiating the determining equation w.r.t. U(o.2) we arrive at 

B{n, m + 1; a, /3) - B{n + 1, m + 1; a, /?) = , 

which gives 

B{n, m; a, /?) = b{m; a, /?) . 

Again, by differentiation, it is shown that the coefficient of a(n; a, (3) and &(rn; a, (3) are independent 
of and U(i respectively. Moreover, in the generic case, the rank of the matrix Q implies 

that the functions ipi are independent on their fourth arguments, i.e. m, U(o,o)7 '^'(i,o)! = 

ipi{n,m,U(^Q Qy, a, P) and 'ip2in,m,U(^0 Q),U(^Q iy, a, P) = 'ip2{n,m,U(^Q Qj;a, P). In this manner, we arrive 
at the following form of the determining equation 

f r m t , ^ hi- \2 a /G("(l,0)j"(0,l)) 

{a{n;a,P) - a(n + l;a,P)) /i(w(o,o), M(i,o)) C„ [— ^ 

' \ "l^i(0,0)j "(1,0)J 



I fhf »\ hi ,1 »\\ h( ^2^) /G(u(i^o),W(o,i)) 

V ^("(0,0), "(0,1)) 

-G(M(i^o),W(o,i))V'("-,"i,U(o^o);a,/3) + /i(u(o.o), '«(o,i))V'("- + 1, "i, U(i^o); /^) 



(32) 



+/i(u(o^o),U(i,o))'0('^,TO + l,M(o,i);a,/3) = + 1, m + 1, a, /?) , 

where we have set 

V'i(n,m,U(o,o);a,/3) +'02(n,rn,U(o,o);a,/3) = ^("-, »77,, U(o,o); «, /?) , 

since they appear additively in the characteristic R. 

The form of the function ip is obtained in the same way as the one used to obtain the general form 
of the characteristic of a Lie point symmetry generator. The substitution of the relevant form of the 
function ijj into the determining equation l|32p and the usage of equation Q = to eliminate U(i i) in the 
resulting equation yield a polynomial in U(q q), q) and ""(q^i). Setting the coefficients of the different 
monomials equal to zero, we come up with an overdetermined linear system of difference equations 
for the unknown functions a{n;a,P), b{'m;a,P) and the functions Ai{n,m; a, P), which occur in the 
general form of the function ip. The general solution of this system delivers the five-point generalized 
symmetries, as well as all the three-point generalized and Lie point symmetries. Thus, equation ([32]) is 
the most general equation for determining the symmetries of a two-dimensional lattice equation, under 
the specific assumptions. The above analysis summarizes to the following 

Proposition 4 Consider the equation (5(m(o,o) j ^(i.o)? ^^(0,1) i ^(1,1)! ck; P) = 0; where function Q is affine 
linear and possesses the D4 symmetry. In the generic case, where the matrix Q defined by 128\) has rank 
3, the characteristic ?ti, ujq q), u^^, 0)1 "(0,1)1 "'^(-i.o); ''^(o, -1)5 Q?; /3) of a five-point symmetry generator 
has necessarily the form 

ip{n,m,uiQ o\;a, P) 

R = a{n; a, /?)P(m(o,o), "(1,0), U(_i^o); P)+b{m; a, /3)P(u(o,o), U(o,i), U(o -1); a)-\ ^ , 

where 

m Hu,x;a,P) 1 

P{u, X, y; a, P) = ^! "> ' 

x-y 2 

and the functions a{n; a, P), b{m]a,P) and ip{n,m,U(^QQy, a, P) satisfy equation 
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7 Extended symmetries on the lattice parameters 

It has been observed [20] that certain lattice equations admit compatible systems of difFerential-difFerence 
equations, with the lattice parameters a, (3 playing the role of the additional (continuous) independent 
variables. Such a compatible system arises from the invariance condition of the lattice equation under 
the action of a generalized symmetry transformation extended to the lattice parameters in a specific 
way. 

As an example, let us consider the lattice potential KdV equation 

("(0,0) - "(1,1)) ("(1,0) - "(0,1)) - a + /3 = . 
Obviously, this equation is not invariant under the scalings 

("(0,0), "(1,0), "(0,1), "(1,1)) (e''"(o,o),e''w(i,o),e''u(o,i),e^W(i_i)) . 
However, if it is further assumed that the parameters a, [3 change according to 

{a, (3) {e^'a,e^' (3) , 

then the equation does remain invariant. In other words, the lattice potential KdV admits the symmetry 
generated by 

"(0,0) 5«(o.o) + 2 ada + 2 [3dfj . 

These observations make it clear that it is useful to extend our considerations to symmetry trans- 
formations acting on the lattice parameters, as well. In this spirit, the present section is devoted to 
symmetry transformations of the equations under study, which are generated by vector fields of the 
form 

w = i?(n, m, U(o,o),"(i,o),"(o,i),"(-i,o),"(o,-i);a,/^)5«(o,o) + ^("■,"^;",/3)^a + C,{n,m;a, (3) dp . (33) 

Acting with the prolonged symmetry generator on the equation and following the analysis in the 
preceding sections, we find, in the generic case, that the component in the w-direction of w takes the 
form 

a(n, m; a, /?)P(u(o,o), "(i,o), "(-i,o); a, (3) + h{m\ a, /5)^'("(o,o) , "(o,i), "(o,-i); ") + ^■0(", m, U(o,o); a, /?), 
where 

"(", X, y; a, [3) = 7;i^.x[u, x\ a, 13) , 

x-y 2 

and the functions a(n;a,/3), 6(m;a,/?), ■(/'("', 'ti, U(o,o); ck, /?), '?('^, ™; ci, /3) and C{n,m]a, (3) satisfy the 
determining equation 

^ I- a\ ^ , 1 «^^ ht \2 n /G(u(i,o),"(o,i))\ 
{a(n;a,(3) - a(n + l;a, (3)) /i(w(o,o), "(i,o)) "uno) T7 ^ 

' ' V ^("(0,0), "(1,0))/ 

I a\ u ,1 Q\\ 1,1 \2 Q /G(u(i o),U(o,i)) 

+ {b(m;a,(3) - b{m + l;a, (3)) /i(w(o,o), "(o,i)) c/«(o n T? ^ 

V "-l"(o,o), "(o,i)J 

+G(u(i^o),"(o,i))V'(n,m,U(o,o);a,/?) + /i("(o,o), "(o,i))'0(" + 1, m, U(i^o); a, /?) 
+^("(o,o),"(i,o))'0(»^,™ + l,"(o,i);a,/3) - Q^„(i.i)V'(" + l,"(i,i);a,^) 
- 2 Q,n(i,i) ^{n, m; a, (3) - 2 Qj3 Q,^^^ C,{n, m;a,f3)=0. 
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Again the form of the function -ip is obtained in the manner described in the previous section, by 
taking into account that 

which foUows from the hnearity of the function Q. The general solution of the above determining 
equation will give us the Lie point symmetries, the three- and the five-point generalized symmetries. 



8 Symmetries of lattice equations with the consistency prop- 
erty 

In this section, we apply the results of previous sections to the symmetry analysis to the class of integrable 
nonlinear equations obtained by Adler, Bobenko and Suris [I] recently. Our results are applicable to 
the ABS equations, because the latter are affine linear and possess the symmetries of the square. They 
are not linearizable and satisfy the condition rankC/ = 3, where G is given by ([28|) . Moreover, these 
equations are integrable in the sense that they satisfy the three-dimensional consistency property. The 
ABS equations are given by entries (Hl-3) and (Ql-4) of the following list. 

(HI) (u(o,o) - "(1,1)) ("(1,0) - ""(0,1)) - a + /3 = (34) 

(H2) (u(o,o) - M(i,i))("(i,o) - ■"(04)) + a)(w(o,o) + W(i,o) + U(oa) + "(i,i)) 

-a^ ^0 (35) 

(H3) a(u(o,o)"(i,o) + '"(o,i)"(i,i)) " /5(w(o,o)W(o,i) + "(i,o)U(ia)) + 5{ct'^ - Z?^) = (36) 

(Ql) a(u(o,o) - W(o,i))(w(i,o) - "(1,1)) - /3(w(o,o) - "(i,o))("(o,i) - "(i,i)) 

+S^aP{a - /3) = (37) 

(Q2) a(u(o,o) - "(04)) (-"(1,0) - ""(i,!)) - /?("(o,o) - "(i,o))("(oa) - U(ia)) + 

a/3{a - /5)(u(o,o) + "(i,o) + W(o,i) + "(la)) - af3{a ~ /3)(a^ ~ af3 + p'^) = (38) 

(Q3) (/3^ - a^)(u(o^o)"(ia) + w(i,o)'^(o,i)) + /'("^ " l)(u(o,o)W(i,o) + U(oa)U(i,i)) 
/«2 ^ <5'(a2-/3')(a2-l)(/32-l) 
-a{(3^ - l)(u(o.o)W(o,i) + "(i,o)U(i,i)) ^ — — = (39) 

(Q4) aoM(o,o)W(i,o)'^(o, 1)^^(1,1) 

+ai(w(o,o)"(i,o)W(oa) + "(1,0)^(0, i)U(ia) + U(o,i)"(i,i)"(o,o) + "(i,i)"(o,o)"(i,o)) 
+a2("(o,o)U(ia) + "(i,o)"(oa)) + a2(M(o,o)"(i,o) + "(oa)W(i,i)) (40) 

+a2(M(0,0)"(0,l) + "(1,0)"(1,1)) + «3(U(0,0) + "(1,0) + "(0,1) + "(1,1)) + 04 = 

Here a, (3 are the lattice parameters and the a^'s in (Q4) are determined by the relations 

ao = a + b , ai — —a/3 — ba , 02 = af]'^ + ba^ , 

ab{a + b) ^ /o 2 92., ~ ab{a + b) 2 /ofl2 92s 

93 92 92 
as = yao - — ai , 04 = —ao - g^ai , 
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with 

= r{a) , = r{(3) , r{x) = Ax^ - g2X - 53 . 

In fact, the results of Sections miH] allow us to obtain edl the Lie point symmetries, the three- and 
five-point generalized symmetries of the ABS equations. The generators of these symmetries are given 
in the following list. We give the corresponding generators, using the symbol of each equation employed 
in the previous list. 

• HI 

Point symmetries : Xi = 9«(o.o) > ^2 = (-l)"~™i9«(o.o) ' ^3 = (-l)"~"'«(o,o)9n(o 0, . 
Three-point generalized symmetries : 

1 U(Q Q\ 

"(1,0) - -"(-i.c) 2(a - P) 

1 U(Q Q) 

■"(0,1) ~ ""(0,-1) 2(a - fi) ' • ' 

• H2 

Point symmetries : Xi = (— . 

Three-point generalized symmetries : 

^^(1,0) + 2u(o^o) + "(-1,0) + 2a _ 2u(o,o) + /9 



Vl 



^(1,0) - "(-1,0) 



-du^o^o) ' V2 = nvi 2(a - (3) ^"<"'°' ' 



^1,0) + 2u(o,o) + "(0,-1) +2/3 _ 2m(o,o) + a 



"(0,1) - "(0,-1) 2(a - P) ' 



• H3 

1. (5 = 0. 



Pomi symmetries : Xi = "(o,o)i9u(o.o) ; ^2 = (-l)""'''""(o,o)(?«(o,o) ■ 
Three-point generalized symmetries : 

_ "(0,0) ("(1,0) +"(-1,0)) a _ "(0,0)("(0,1) +"(0,-1)) a 

Vl — ^uro 0) ' '^2 — C'tijo 0) • 

"(1,0) - "(-1,0) "(0,1) - "(0,-1) 



Five-point generalized symmetries : w = n vi + m V3 . 
2. 6^0. 

Point symmetries : xi — (— 1)"+'"m(o,o)5m(o q, . 

Three-point generalized symmetries : 

_ "(o,o)("(i,o) +"(-1,0)) + 2a<5 _ "(o,o)("(o,i) +"(o,-i)) + 2/3<5 

Vl — C«{0,0) ' V2 - C„( 

"(1,0) - "(-1,0) "(0,1) - "(0,-1) 



Five-point generalized symmetries : w = nvi -I- mv2 — '2'°'' ^"(o o) 
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• Ql 



1. 6 = 

Point symmetries : xi = M(o,o)^«(o,o) ' ^2 = '"(o,o)5«(o,o) ' ^3 = . 
Three-point generalized symmetries : 

(U(i,o) - ■«(o,o))(w(o,o) - W(-1,0)) „ 

■^1 = '^"(0,0) ' 

"(1,0) - "(-1,0) ^ ' 

("(0,1) - "(0,0)) ("(0,0) - "(0,-1)) r. 

•"(O.l) - "(0,-1) ' ' 



Five-point generalized symmetries : w = n Vi + m V2 . 

2. ^7^0 

Point symmetries : xi = du(o,o) ■ 

Three-point generalized symmetries : 



2r2 



("(1,0) - "( o,o))("(o,o) - "(- 1,0)) + a (5 

"(1,0) - "(-1,0) 



*1 ~ '^"(O.O) 



("(0,1) - "(o,o))("(o,o) ~ "(0,-1)) + /3 ^ Q 

"(0,1) - "(0,-1) ' ' 



Five-point generalized symmetries : w = nvi + mv2 — "(0,0) c^m(o,o) • 

• Q2 

Three-point generalized symmetries : 

("(1,0) - "(o,o))("(o,o) - "(-1,0)) + ("(1,0) + 2u(o,o) + "(-i,o))a^ ~ o 

"(1,0) - "(-1,0) ' ' 

("(0,1) - "(0,0))("(0,0) - "(0,-1)) + ("(0,1) + 2u(o,o) + "(0,-1))/?^ ~ a 

"(0,1) -"(0,-1) 



Five-point generalized symmetries : w = nvi + mv2 — 2u(q g) . 
• Q3 

Point symmetries : If ^ = 0, then it admits one point symmetry with generator xi = U(o,o)9i 
Otherwise, there are no point symmetries. 
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Three-point generalized symmetries : 

_ 2Q(a2 + l)u(o^o)(w(i,o) + W(-i,o)) - 4a2(u2^ + U(i_o)U(-i,o)) - ("^ - "^^^'^^ « 

W(1,0) - ■"(-1,0) 

_ 2/3(/32 + l)u(o,o)(^^(o,i) + ^i(o,-i)) - 4/3^("(o,o) + ^(o,i)W(o,-i)) - (Z?^ - l)^'^^ 

^2 — ■ "tiro.o) ■ 

M(0,1) - U(0 -1) 



• Q4 

Three-point generalized symmetries : 

_ ("(1,0) - '"(-i,o))/.«(i,o) (^(0,0). "(1,0)^ - 2/(-»(o,o), -»(i,o), g) ^ 

W(1,0) - "(-1,0) 

_ ("(0,1) - ■«(0,-l))/,«(o,i)("(0,0),"(0,l),/?) - 2/(m(o,0),"(0,1),/?) 

"(0,1) - "(0,-1) 

where 

/("(o,o),"(i,o),a) = (^"(o,o)"(i,o) + a(u(o,o) +"(1,0)) + -("(o,o)+"(i,o)+a)(4au(o,o)"(i,o)-53) • 

The extended symmetry transformations acting on the lattice parameters along with the corre- 
sponding determining equation have been presented in Section [71 Using these results, we find that the 
integrable lattice equations of Adler, Bobenko and Suris admit the extended symmetries of the following 
list. 

• HI 

Point symmetries : X4 — U(^Q Q-)d.a^f, q, + 2ada + 2/35^5, X5 = da + dp . 

Three-point generalized symmetries : 

V5 = A{n)^^l + {A{n) - A{n + 1)) vg = S(m)v3 + (S(m) - B{m + 1)) dp. 

• H2 

Point symmetries : X2 = "(o,o)f^U(o.o) + "f^a + Pdp, X3 = dufg a^ - 2da - Idp . 

Three-point generalized symmetries : 

V5 = A(7i)vi + {A{n) - A{n + 1)) vg = S(m)v3 + {B{m) - B{m + 1)) dp. 

• H3 

1. (5 = 0. 

Point symmetries : X3 = q;(9q + /3i9/5 . 

Three-point generalized symmetries : 

V3 = A{n)vi - (^(n) - A{n + 1)) ctc)„, V4 = B{m)v2 - {B{m) - B(m + 1)) pdp. 
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2. 5 7^0. 

Point symmetries : X2 = 'W(o,o)5«(o,o) + 2q!9c« + 2/39/5 • 

Three-point generalized symmetries : 

V3 = ^(n)vi - {A{n) - A{n + 1)) ad^, V4 = B(m)v2 - {B{m) - B{m + 1)) 

• Ql 

1. (5 = 0. 

Pomf symmetries : X4 = + /39/3 . 

Three-point generalized symmetries : 

V3 = A(n)vi - (A(n) - A{n + 1)) aa„, V4 = B(m)v2 - (B(m) - B{m + 1)) 

2. 5^Q. 

Point symmetries : X2 = 'U(o,o)^«(o,o) + + f^dfj. 
Three-point generalized symmetries : 

V3 = ^(n)vi - (^(n) - A{n + 1)) aQ^, V4 = B{m)v2 - {B{m) - B{m + 1)) /35/3. 

• Q2 

Point symmetries : xi = 2u(o,o)<?u(o 0) + '^^a + /5^/3 • 

Three-point generalized symmetries : 

V3 = >l(n)vi - (A(n) - A(n + 1)) ada, V4 = B{m)v2 - {B{m) - B{m + 1)) pdp. 

• Q3 

Three-point generalized symmetries : 

V3 = A(n)vi - 2 (^(n) - A{n + 1)) a2(Q,2 _ 1)5^ ^ 
V4 = B{m)v2 - 2 (B(m) - B{m + 1)) /32(/32 - 1)00 . 

• Q4 

Three-point generalized symmetries : 

V3 = A(n)vi + (A(n) - A(n + 1)) {4a^ - - 53) , 
V4 = B(m)v2 + (B(m) - B{m + 1)) (4/3^ - g2p - 53) 9/? . 

9 Symmetry reductions to discrete analogues of the Painleve 
equations 

In the same way as in the case of the partial differential equations, the symmetries of a partial difference 

equation provide an effective means of constructing (whole classes of) special solutions. They are the 
ones that retain their form when acted on by some of the transformations leaving the equation invariant. 
Hence, they are referred to as group invariant solutions. 
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Group invariant solutions of certain integrable lattice equations are known to be associated with 
integrable mappings and discrete versions of the Painleve equations. Reductions to integrable mappings 
first appeared in [22j . where the periodic boundary value problems for the lattice potential KdV were 
studied (see also [6l [25]). On the other hand, reductions to discrete Painleve equations originally 
appeared in [20j . where the higher symmetries of the latter equation were exploited. For a more recent 
account on the subject we refer to [101 1121 flQl l2T]. 

Here, we use the results obtained above to study group invariant solutions of the class of lattice 
equations introduced in Section |31 as well as initial value problems associated with them. More specifi- 
cally, we relate the existence of symmetry reductions to specific Cauchy problems, leading to a unique 
solution. The existence of such group invariant solutions essentially follows from the fact that every 
equation in the class admits a five-point generalized symmetry with generator v„ -I- Avm (see Proposition 
[3]). In the generic case, these reductions lead to four dimensional mappings. We show that the existence 
of a Lie point symmetry of the lattice equation, which is also compatible with the five-point symmetry 
constraint, suffices to reduce the resulting mapping to a third dimensional one. As an illustration we 
apply the results just mentioned to the lattice potential KdV (HI). 

9.1 Symmetry reductions and Cauchy problems 

Let Q = be a partial difference equation and v = R[u\ du^„ an infinitesimal generator of a symmetry 
transformation acting on the space of the dependent variable u, only. The square brackets in the 
symmetry characteristic R[u\ means that, the latter is, in general, a function of (n, to, M(o.o)) and the 
shifted values of u in both directions of the lattice up to some order fc, i.e. R : 1? x ^ C for 

some fixed fc e IN. 

Definition 1 A function u : 1? i-^ G is called an invariant solution of the lattice equation Q = 
under the symmetry v, if it satisfies the lattice equation Q = and the compatible constraint v(u) = 
(equivalently R[u\ — 0). 

A natural question is whether non-trivial solutions of this kind do really exist. Consider a generic 
lattice equation Q = of the form ([6]) , where Q satisfies the properties of Section [S] According to 
Proposition |3l the lattice equation Q = always admits two generalized symmetry generators v„, v^. 
Thus, the invariant solutions under the symmetry generator 

satisfy the lattice equation and the compatible symmetry constraint 

C[u] := i?(w(o,o),'«(i,o),W(-i,o);a,^) + A i?(u(o,o), "(o,i)' "(o -i); ") = 0: (41) 

where 

R{u,x,y;a,/3) = ^(^''^'"'^) _ x; a, /5) . 

x-y 2 

Let us now consider the following initial value problem. Given the values of u assigned on the four 
points depicted with black in Figure 31 one can determine uniquely the values of u along three vertical 
columns of the square grid, by using the lattice equation (O points) and the symmetry constraint (-1- 
points). It should be noted that, due to the affine linearity of the function Q, the symmetry constraint 
(|41l can be solved uniquely for each one of the values q), q), U(o,_i) and U(^Q iy Then, the values 
of u on the remaining vertices of the lattice can be determined, by using successively the symmetry 
constraint. In this setting there will be points where the values of u can be determined by using either 
the lattice equation, or the symmetry constraint (® points). However, since the symmetry constraint 
is by definition compatible with the lattice equation, the values determined by the two different ways 
described above are, necessarily, identical. Thus, for generic initial data, such a group-invariant solution 
exists and is unique (see the discussion in PIUUH^I])- Let it be noted however, that such a symmetry 
reduction leads to four-dimensional mappings. 
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Figure 4: A Cauchy problem on the lattice. Figure 5: The Cauchy problem on the dual (chess- 
board) lattice. 

Nevertheless, under certain conditions, symmetry reductions of Cauchy problem of the type under 
consideration lead to mappings which are three-dimensional. To see how this is obtained, let us require 
that a given Lie point symmetry of the lattice equation, is also a symmetry of the constraint C[u\ = 0. 
If X = X(n, m, it(o,o)) '^u^o o) is the generator of the Lie point symmetry, then the above requirement is 
equivalent to the following commutation relatior0 

[x,vc]=0, mod C[u] = 0. (42) 

The existence of the point symmetry generated by x allows one to define lattice invariants {reduced 
variables) assigned on the edges of the square grid, in terms of which, the lattice equation can be casted 
as a map. The dynamics of this map is restricted to the white squares of the chessboard in Figure [5l 
Moreover, since x is also a symmetry of the constraint C[u] = (defined now on the vertices of the 
shaded squares of the chessboard), the latter can be written in terms of the invariants of the vector 
field X on J^^'"^^ In this way, the corresponding invariant solutions are constructed from the associated 
Cauchy problem with initial values xi,X2,X3 assigned on the edges (Figure [5]), and, consequently, the 
mappings obtained by this procedure are, in general, three dimensional. 

A specific example illustrating the case of reductions to three-dimensional mappings is provided by 
the lattice potential KdV equation (l34|) . The most general symmetry constraint on a cross configuration 
of points is obtained from a linear combination of the symmetry generators presented in the previous 
section, i.e. 

3 4 

vc = ^ Xi + ^ Vi , (43) 

where Xi, Hi are arbitrary complex parameters. Here, {xi,X2,X3} is the set of Lie point symmetry 
generators, which span a Lie algebra q isomorphic to 5o(l, 1), and {vi, V2, V3, V4} is the set of the three- 
point generalized symmetries . We note that, under the transformation u 1— > (— 1)"'+'"m, the symmetry 
generator X2 is mapped to xi. Thus, for the symmetry reduction using the invariants of one-dimensional 
subalgebras of g, it is sufficient to consider two inequivalent cases, namely x = xi or X3. 

9.2 Symmetry reduction of HI using the invariants of Xi 

We consider first the case x = xi and recast the lattice equation HI as an invertiblc map TZ : CP^ x 
CP^ f-> CP^ X CP^ To this end, we introduce the following lattice invariants along the orbits of xi 

V(o,o) = W(i,o) - W(o,o), ■"'(0,0) = W(ia) - "(1,0), ■'^(o,!) = "(14) ~ "(0,1)' "^(-i.o) = "(0,1) ~ "(0,0), (44) 

^Indeed, the commutation relation [x,V£;] = 0, yields x(C[«]) — C[ii]9t,jp {X(n, m, U(q qj )) = which, in virtue of 
C[u] = 0, means that x is also a symmetry generator of the constraint C[u] = 0. 



20 



where w, 



G jI'^^' and J*"^^ is the space obtained from J^^^ by removing (u(2_o)j 



t(2) 



There is a functional relation among the above invariants, namely 



(2,0)j "(0,2)J 



(45) 



following from the fact that the space of invariants along the orbits of x on is three dimensional. 
On the other hand, since x is a symmetry generator of equation HI, the latter can be written in terms 
of the invariants (f44|) in the form 



(w(o,o) + u'(o,o))(i'(o,o) - W(-ifi)) + r = , 



(46) 



where r = a — p. Equations (|45p . (|46p can be uniquely solved for ^(o,!): w;(-i,o) 
and conversely, implying the existence of the following invertible map 



in terms of v 



(0,0) J ^{0,0) y 



h + h 



/1 + /2 



(47) 
(48) 



where 

(/l,/2,/3,/4) = (W(0,0), W(o,o),W(o,l),'U;(_i^O)) • 

Moreover, equations (P5|) . (|46p can be MnigweZi/ solved for z;(o.i) , u'(o,o) in terms of W(o, 0)1 if(-i.o)7 implying 
also the existence of the invertible map 



^(/l,/4) = (/3,/2) 



-h 



h-h 



(49) 



The existence of the above invertible map 1Z is associated with the quadrarationality property of the 
birational map TZ. Maps with this property were studied recently in [5], in connection with the Yang- 
Baxter relation. Here, the quadrarationality property of the map R allows one to uniquely determine 
two of the four values fi assigned on the vertices of the white squares of the chessboard in terms of the 
remaining two, as shown in Figure [6l 

fs /a 



/2 f4 



h fl 

Figure 6: The invertible maps TZ, TZ. Their dynamics is restricted to the white squares of the chessboard. 

We proceed by writing the symmetry constraint vc(u) = 0, where vc is given by equation (|43p . in 
terms of the invariants (UH) and their shifts. To this end, it is necessary that relation ([1^ holds. A 
direct calculation shows that 

[xi,vc]=A3X2 + i^^^^Xi, (50) 
2 a — p 

from which we conclude that the parameters Ai, /ii should be rectricted by the following relations 

A3 = , ^2 = • (51) 

Taking into account the above relations, one can write the symmetry constraint vc(u) = in terms of 
the invariants (|44p in the form 





C(n,m,5ii,52,53,54) := Ai + A2 (-1) ^ H H ■ = 0. 



gi + 92 



93+04 



(52) 
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where 

{91,92,93,94) = {v(ofi),V(-iM),W(-i,o),W(-i^^i)) . (53) 

We consider now the Cauchy problem on the chessboard with initial values (xi, X2, X3), as shown in 
Figure [5l The updated values {x'i,X2,x'z) in the n-direction of the lattice can be found by using the 
invertible maps TZ, TZ and the symmetry constraint ([52|l . Indeed, the value X2 is the second component 
of the map 'R{x3,X2), and obviously Xg — xi. The updated value X3 is computed in three steps. First, 
we find X4 by solving the equation C(n, m, X3, xi, X4, X2) = for X4, next from the map i?~^(x3,X4) = 
(x5,X6) we pick the second component xg, and finally the value Xg is found by solving the equation 
C{n + 1, m, X3, X3, xg, X2) = 0, for the corresponding variable. A straightforward calculation shows 
that the updated values (x'l, X2, x'^) are found from the following non-autonomous system of first order 
ordinary difference equation in the variable n, 

x'l = X3, (54) 

T 

x'2 = -X3 H , (55) 

X3 - X2 

2-' _ rcin + l,fii) 

X3 — X3 + ,^ \ ' ' 

( X1+X3 + + ) (^2 - X3)2 + c(m, fis) (x2 - X3) - rX{n + to + 1) 

where c{n, ^) = ii2'n + ji and A(n) = Ai + A2 (—1)". The above system can be decoupled for the variable 
y{n) := X2 + X3 = X2(ri + 1, m) + X3(n, to) leading to the second order difference equation 

— c(n+l,^i) + c(TO,//3) + y(n) A(n + TO.+ 1) — A(n + m), (57) 



j/(n + 1) y(n) + r y{n)y{n — 1) + r ' ' ?/(n) 

which is known as the asymmetric, alternate discrete Painleve II equation [S] , [llj . 

9.3 Symmetry reduction of HI using the invariants of X3 

We now turn to the case where we use the invariants of x = X3. Following the preceding considerations 
(see also [H]), we introduce the following invariants on jI^-* along the orbits of X3 

^'(0,0) = ""(1,0) U(o,o), "("(CO) = "(1.1) "(1,0), "^(0,1) = "(1,1) "(0,1), "'(-i.o) = "(0,1) "(0,0) • (58) 
The above invariants are functionally related by 

"(0,0) "(0,1) = ^(0,0) "'(c-i) • (59) 
On the other hand, equation HI can be written in terms of the invariants (j58p in the form 

"(0,0) - "^(-1,0) - "'(CO) + "(0,0) = r . (60) 

Using equations ([CT]) . ((M)) we obtain the invertible maps TZ : (/i,/2) ^ (/sjA), T^[fi, fi) ^ (/3,/2)- 
Explicitly they are defined by 



(61a) 



7^(/l,/2) = (/3,/4)= f./2(l + ^r^),/i(l + ^^)V 

V /2 - /i h- h J 

n{h,h) = nfi,h) - (/3,/2), (61b) 

where {fi,f2,fz,U) = ("(0,0) , ""(0,0), "(0,1) , "'(-1,0))- Next we write the symmetry constraint vc(u) = 
in terms of the invariants ([55]) . where vc is given by A straightforward calculation leads to the 

following result 

[x3,vc] = -A1X2 - A2X1 . (62) 
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Hence, relation (j42|) holds whenever 

Ai^O, A2 = 0. (63) 
In this setting, the symmetry constraint vc (u) = is written in terms of the invariants ([58|) as follows 

C{n,m,gi,g2,g3,9i) ■■= Xsi-l) H ;^ h \ = 0, 64) 

2r 51 - 52 53 - 54 

where 

(51,52,53,54) = (i'(o,o),i^(-i,o),W(_i,o),W(_i _i)) . (65) 

The corresponding Cauchy problem on the chessboard is reduced now to the solution of the ordinary 
difference equations defined by the mapping M : {xi,X2,X3) i-> {x[,X2,x'^). The updated values are 
found in the similar manner as in the previous case where now the invertible maps TZ, TZ and the 
symmetry constraint are given by equations (pT|) . (|64p . respectively. More precisely we have 

x'l = X3 , (66) 



4 = X3 1 + , (67) 

X2-X3J 



c(n + 1) ^ + c(n) = c(n + 1) + d{m) + X{n + m){x2 - X3) + X{n + m + l)(x2 - X3), (68) 

X3 — X\ X3 — Xl 



and X3 is found by solving 

C3 - X2 , . X3 - X2 

where 



c{n) = ^2n + fii , rf(m) = ^4TO + ^3, A(n) = A3(-l)" + . (69) 

2r 



If /i2 = /i4 then equation (j68p can be integrated once. Due to a compatibility condition, the arbitrary 
function of m in the discrete integration is specified up to a constant, yielding the following result 

c(n) ^^-^ = (-1)"+™ (A3 X2+p) + i(c(n) + d{m)) + ^ , (70) 

where p is the complex constant of integration. One may solve equation (j70p for X2 and use the result 
to eliminate X2 from equation (|67p . using (p^. This leads to a second order difference equation for x\, 
which we omit because of its length. 

10 Conclusions and perspectives 

We have presented a symmetry analysis of a class of lattice equations on 1? ^ which are characterized by 
affine linearity and D4-symmetry. Once a specific equation in the class is given, the results summarized 
in Propositions [3] and [4] explicitly determine the characteristics of its three- and five-point generalized 
symmetry generators. Applied to the integrable equations obtained in classification 1^ these results 
allowed us to determine all Lie point, three- and five-point generalized symmetries admitted by the 
above equations. The results obtained in this fashion provide a proof that these lists of symmetries 
are exhaustive. From this point of view, the present work constitutes a generalization of the studies on 
these equations presented in [15], [27], [28], [30] . 

The effectiveness of using Lie point symmetries of integrable lattice equations in obtaining Yang- 
Baxter maps was demonstrated in [53], [53]. In particular, it was shown there that the Yang-Baxter 
variables can be chosen as invariants of the multi-parameter symmetry groups of the integrable lattice 
equations. Here, it was shown that this connection, combined with the quadrirationality property of the 
associated Yang-Baxter maps [2j, can be used in obtaining group invariant solutions of a lattice equation. 
Specifically, we considered a Cauchy problem for the lattice potential KdV equation compatible with 
the most general symmetry constraint on five points. The solution of this initial value problem was 
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constructed by solving a second order mapping, which represents a discrete analogue of the Painleve 
equations. 

The preceding symmetry analysis may also be applied to those equations in the class considered 
above which do not have the consistency property. It would be interesting to consider the corresponding 
symmetry reductions and investigate various properties of the resulting mappings, such as the singularity 
confinement [S^ and the algebraic entropy [4j. Moreover, it would be interesting to extend the results 
obtained here to lattice equations on 1? ^ which possess symmetries other than the symmetries of the 
square. Work in this direction is in progress. 
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Appendix. Proof of Proposition [3] 

Let 

1 1 

v(2) = ^^i?(n + i,r7i + j,U(,j),M(j+i_j),U(,_i_j);a,/3)9„,.^.) 

be the second prolongation of the generator v of a three-point generalized symmetry. In this case, the 
infinitesimal symmetry criterion takes the form 

<9,U(o,i)-R(", m + 1, ^(o,!), a, /?) + Q.^d,!, + 1, m + 1, ^(24), U(o4); ^) 

+Q,ti(o,o)-^('^>™>'"(o,o),'«(i,o), W(-i,o);a,/3) + Q,«(i,o)-R("- + l,"^:"(i,o),M(2,o),"(o,o);a:/3) = 0. (71) 

The last equation must hold on every solution of the equation 

Q = (3(w(o,o),"(i.o),'"(o,i),U(i,i);a,/3) = 0, 

and its shifted consequences, i.e. 

g = (9(M(-1,0),U(0,0), "(-14)' "(0,1)5 = 0, Q = Q(M(1,0),M(2,0),'"(1,1)''"(2,1);"i/5) = 0. 

Using the above equations, one may express the values ^(24), U(o,o) and U(_i.o) in terms of the remaining 
ones. 

One can eliminate the shifts of the characteristic R from equation (1711) by differentiating the latter 
w.r.t. U(_i 1), which yields 

^.^(o,!) -^,"(-1,1) ™ + 1, ""(0,1), ^(-i,!); a, /?) 

dui_i Q\ 

+ Q,u,o 0) ^,"(0 -1) ("> "(0.0), U(i.o), "(-1.0); a, P) J. — = , (72) 

C"(-l,l) 

since only U(_i.o) implicitly depends on 1) through the equation (5 = 0. Next, we divide equation 
([72)1 by and take the total derivative of the resulting equation w.r.t. it(i,o), i-e. 

D«(i.o) ( (n,m,U(o,o), "(1.0), "(-1,0); «,/3) ] =0, (73) 

V^>"(o,i) C'"(-l,l) / 
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■"(-1,0) "(0,0) "(1,0) 

Figure 7: The points of the lattice and the corresponding equations. 



where 



du, 



du 



(0,0) 



du 



(1,0) 



du 



(0.0) + 



5u(_i,o) 5"(o,o) 



(0,0) 9zi(i,o) ■ 



Writing equation (|73p exphcitly, one arrives at 



(G',«(i,o)("(i,o),w(o,i)) + (M(_1,0),'«(0,1))) R ,«(-!, 0) + G(w(_i,o),"(o,i))^ 

'*,"(i,o)('"(0,0),'"(l,0)) 



+ G'(w(i,o),U(o,i)) R 



^("(0,0), ■"(1,0)) 



(74) 



^("(0,0), "(0,1)) ^, 



'."(O.OjMC-l.O) 



"(0.0) ("(0.0)' "(1,0)) 
/i("(0,0),"(l,0)) 



-^,"(-1.1 



0, 



where we have omitted the arguments of the function m, U(o_o); "(i.o); "(-1,0)5 Q^i /^) for simphcity, 
and the equations Q = 0, Q = have been taken into account to evaluate the derivatives of U(_2,o) ^-nd 
"(0,0): i-e- 



3m(0,0) 


_ ^."(i.n) _ 


^("(0,0), "(0.1)) 




G(u(i,o),U(o,i)) 


9W(1,0) 


*3,"(o,o) 


G(w(i_o),U(o,i)) ' 


9u(o,l) Q,"(o,o) 


^("(0,0), "(1,0)) ' 


9"(-l,0) 


§."(-1.1) 


G(-u(_i^o),"(o,i)) 


f^"(-l,0) §,"(0,0) 


G(u(-i,o),"(o,i)) 


^"(-1,1) 


§■"(-1.0) 


/i("(-l,l),"(0,l)) 


5"(0,0) §,«(_!, 0) 


^("(0,0), "(0,1)) 



Next, we substitute the derivatives of the polynomials G appearing in Equation ([7^ by the relation 

G'(^(1,0);"(0,1)) - G'(m(_i^o),M(o,i)) 



^,"(1.0) ("(1,0): "(0,1)) + G,.U(_i_o)(u(_i,o),"(o,i)) = 2' 



"(1,0) 



^(-1,0) 



which follows from the fact that the polynomial G is quadratic and symmetric in its arguments. Upon 
these substitutions, equation ([7^ simplifies to 



- 'l("(0,0))"(0,l)) ^^,"(o.o)«(-i. 



+ ^(■"(1^0)7 "(0.1)) ^-R.n(i,o)«(-i,o) ~ ~ 



."(0.0) ("(0.0), "(1,0)) 
°' ^("(0,0), "(1,0)) 

^,"(i,o)("(0,0)j"(l,0)) 



R 



"(-1.0) 



(-1,0) 



i(O.O) 



+G'(-U(_i_o),"(o.i)) ^^,«(_i,o)"(-i,o 



^'("(0,0), "(1,0)) 

2 



^,"(-1.0 



(75) 



i(l,0) 



^(-1,0) 



"^,"(-1,0) I — • 



The last equation involves the values of u assigned on four vertices of the lattice (the black ones in Figure 
[7]) and should hold on every solution of equation Q = 0. Equation ([75)1 depends on M(o.i) through the 
polynomials /i("(o,o)7 "(o,i)), ^(^(1 g), M(o_i)) and G(u(_i q) , "(o,i)), which are in general quadratic in 
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M(o,i). Thus, it is necessary to set the coefficients of different powers of W(o,i) in the latter equation 
equal to zero, i.e. 



'l(u(0,O),0) G'(m(i^o),0) G(u(_i^o)70) \ / ^,U(o.om-i.O) + -^i ^,-u.i-i,o) 



/i'(m(o,o),0) G'(w(i,o),0) G'(u(_i.o),0) 
/i"(«(o,o),0) G"(u(i,o),0) G"(w(_i,o),0) 




A2 R ; 

2R. 



•,«(_1,0)M(-1,0) «(1,0)-M(-1,0) 

where the prime denotes differentiation w.r.t. U(o,i) and 

. ^,«(oo)("(0,0),"(l,0)) . 2 ft,n(i o)("(0,0),'"(l,0)) 

-4i — —, r , A2 









(•) 







(76) 



^("(0,0), ""(i^o)) ' -"(-1,0) - U(0,0) ^(U(0,0),"(1,0)) 

In the generic case where 

rank^ = 3, (77) 

with 



h{x,y) G{x,z) G{x,w) 
G = \ h^,j:{x,y) G^x{x,z) G.^{x,w) 
h^xxi^^y) G^xx{x,z) G,xx{x,w) 



(78) 



x = 



then system (j76p has the unique solution 



n _ ^,"(0.0)('»(0,0)-"(1.0)) n 

-"^."(O.ofc-i.o) — . / N, ^,"(-1,0) ' 

'11^(0,0), U(i,o)j 

P - 2 , ^"(i.o)(mo)^"(i,o)) A „ ,7q^ 

-n.,u(_i.o)tt(-i.o) ^:"(-i,o) ■ 

"(1,0) ^ "(-1,0) 

Let it be noted that, even though the rank condition is violated this is a solution of system (j76p and, 
consequently, of equation (|75p . 

Integrating system ([7^ , we find that the characteristic reads the form 

i?(n,m,U(o,o),U(i,o),'"(-i,o);a,/^) = A{n,m;a,(3) ^("("-O)' "(^-O)^ + r(n, m, M(o,o), U(i,o); /?) , (80) 

"(1,0) - ""(-1,0) 

where A and r are arbitrary functions of their arguments. 

Next, we substitute (|80l) into the determining equation ((7T|) . Using the relatioijl 

L / ^ _L ^ ^ fe(u(o,o),M(i,o)) - ^(m(-i,o),M(o,o)) 
^,"(io)(«(o,o),W(i,o)) + 1 o)("(-i.o),'^(o,o)) = 2 , (81) 

""(1,0) ~ "(-1,0) 

and its shifted versions, we find that it is convenient to set 

, -^("-, w;a,^)/i,„ (m(o,o),U(i,o)) + 0('^,™,M(o.o),W(i.o);a,^) 
r(n,m,U(o^o),'"(i,o);a,P) = ^ ■ 

Upon these substitutions the characteristic takes the form 

m /^ fe("(o,o),'»(i,o)) 1, ^ \\ , ^Af m ^ao\ 

A(n,m;a,/i) x^.^d 01 (^(0,0), "(1.0)) + ™, W(o,o) , U(i,o); a, /?) (82) 

\"(i,o) ^ ""(-1,0) ^ / ^ 



^To prove this relation, we write /i(u(o,o) > -"(LO)) = P2("{o,0))«(i.o)+Pl("(0,0))"(i,0)+Po(«{o,o)) and ?t(M(_i,o) , M{o,o)) = 
<S^ ('^('"(0,0) 1 '^(l,0)))i since the latter function is the backward shift of the former in the n-direction. Using the fact that 
they are symmetric, a straightforward calculation implies the result. 
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and the determining equation becomes 



2/i(u(o,o),U(i.o)) 



/ / 2h{ 

Q,n,oo) A{n,m;a,P) -/i.^n oi ("(0,0), "(1,0)) H 

\ \ ■ "(1,0) ^ ""(-1,0) 

+Q,«(i,o) I + 1' P) I -^:«(2,o) ("(1,0)> "(2,0)) + 



"(2,0) - "(0,0) 



+Q,«(oi) fA(n,m + ^ (u(o,i), "(1.1)) H (^fCHl^JfiMll^ 



(83) 



+ 1,771+ l;a,[3) ^-/i,«(2.i)(w(i,i),M(2,i)) + 



2/l(W(i,i),U(2,l)) 



^(2,1) 



i(0,l) 



- 5„(5,„(0)) 



where we have omitted the arguments of the function (p and its shifted values. Equation (|83p involves 
the values of the function u assigned on the eight vertices of Figure [T] Using the equations Q = on 
the three faces, we eliminate three of these values, and we have chosen to eliminate U(-i,o), 1) and 

"(2,1)- 

The only terms in (j83p which depend on M(_i.o) and M(_i.i) appear in the coefhcients of A{n^ m; a, [3) 
and ^(77, 771 + 1; a, /?), respectively. We take the total derivative of the determining equation (|83p w.r.t. 
"(-i.i)j since q) depends on through the equation (5 = 0. After a lengthy calculation, this 

simplifies to 



{A{n,m;a,l3) - A{n,m + 1;q!,/3)) /i(u(o,o), "(1,0)) ^("(0,0), "(0,1)) G(?7(i,o) , U(o,i)) 



0, 



(84) 



where we have used equations Q = 0, Q = to eliminate 7t(i,i) and 7t(_i 0), respectively, and relation 
(jl4bp . Equation ([M]) implies that A[n,m]a, (3) does not depend on 771, i.e. 



^(77, 7m; a, (3) = 0(77; a, (3) . 



(85) 



Thus, taking into account the latter, the terms in (I83p involving ^(77, 777; a, (3) and ^(77, 777 + 1; a, f3) are 
independent of 7i(-i,i), and consequently they can be evaluated at 7i(-i,i) = 0, simplifying to 



a(n; a, (3) 



/ ( 

Q,Ho,o) ( ~^,"(i,0) ("(0,0), "(1,0)) + 2/l(7i(o,0), "(1,0))- 



,"(1,0) ("(0,0), "(1,0): "(0,1): 0) 
Q("(0,0),"(1,0),"(0,1),0) 

+y,u(o,i) 1 "-,«(i.i)("(o,i),"(i,i)) 

\ "(1,1) 



(86) 



Moreover, it turns out that the determining equation, apart from the function 0, does not depend on 
the value 77(2, o)- Indeed, substituting ([85]) in equation ()83p . the total derivative of the resulting equation 
w.r.t. 77(2,0) is identically zero. Thus, the relevant terms in ()83p . except the function (p, arc independent 
of W(2,o), and hence we can set U(2,o) = in them. These can be written as 



"'*,"(2,0)("(1,0):"(2,0)) + 
-^,«(2,l)("(l,l),"(2,l)) + 



2/7(7^(1,0), 77(2,0)) 
"(2,0) - "(0,0) 

2/7(7^(1,1), 77(2,1)) 
"(2,1) - "(0,1) 



f (2,0)=0 



f (2,0)=0 



(0,0) 



^"(0.0) 



^("(0,0), "(1,0)) 



(87a) 



^0,0) 



'*,"(o,i)("(0,l)>"(l,l)) 



-2/7(77(0,1), 77(1,1)) 



'3,"(o.i) 



Q 



(87b) 
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Substituting ([85| . (|86|) and ([87|) in equation (|83|) . the latter simplifies to 

, ^."(i.o) (W(0,0), U(l.O), "(04), 0) 



a{n; a, (3) 



,"(0.0) ^ -^,«(i,o)("(0,0),'"(l,0)) + 2/l(u(o,o),U(l,0)) 

+'9,"(o. 



'9("(0,0):"(1.0)>"(04)'0) 

2/i(w(o,i),W(i,i)) , , , 
'i,u(i,i)(,"(oa),'"(i,i)J 



+a(n + 1; a, (3) 



"(1,1) 

W(0,0)^,U(o,o) (W(O.O), "(1,0)) - 2ft,(u(o_o),"(l,0)) 



,"(1.0) 



+Q,"(i.i) ( ^,"(0.1) ("(04), "(1,1)) - 2ft.(u(o,i), "(1,1)) 



"(0,0) 

Q,"(o,i)(0, "(1,0), "(0,1), "(1,1)) 



Q(0, "(1,0), "(0,1), "(1,1)) 
+<3,"(o.o)'^('^'™'"(o,o),"(i,o);a,/3) + Q,n(i.o)'/'('^ + l,™,"(i,o),"(2,o);a,/3) 

+Q,«(o,i)'^("-, "1 + 1, U(o,i), "(1,1); a, ^) + + 1, m + 1, U(i4), ii(2,i); a, /?) = . 

It remains to simplify the coefficients of a(n; a, /3) and a{n + 1; a, /?) by eliminating the value "(i,i)- 
For this purpose, we use the relations ()14cp . (ll4dp and the following ones 



"(1,1) (fc,n(i.i)("(0,l),"(l,l)) +G',M(i_o)(M(1,0),U(0,1))) " 2/j.(M(o, i) , U(i, i) ) _ 

2"(i,i)G(u(i,o),"(o,i)) I 

"(0,0) ('i,U(o,o)("(0,0),"(l,0)) + G,„((,_jj(M(i_o),"(0,1))) - 2/l(lt(o,o),"(l,0)) 



'(1.0) 



(0.1) 



2"(o,o)G(u(i,o),U(o,i)) 



«(i.i>=0 



«(o.o)=0 



^^("(0,0), "(l,0))G',,i(i_o)("(l,o),"(0,l)) - ^,"(1,0) ("(0,0),"(1,0))G(m(i,o), "(0,1)) _ Q«(i,i 

^("(0,1), "(l,l))G,U((,_i)("(i,o),"(0,l)) - ^,«(o,i) ("(0,1),"(1,1))G'(W(1,0), "(0,1)) Qm(o,o 

which hold in view of the equation Q — Q. 
Finally, equation ([55)1 simplifies to 



(89) 
(90) 
(91) 



(a(n;a,/3) - a(n + l;a,/?)) /i(m(o.o), "(i,o))^ 5«(i o) fT/^^^'"^ ' "'"'^l'' 

V ^("(0,0), "(1,0)) 

+G'("(i,o),"(o,i))'/'(","^,"(o,o),"(i,o);a,/3) + ''-("(0,0), "(o,i) )'/'(" + 1, "t-, "(i,o), "(2,0); /^) (92) 

+/i(m(o,o) , M(i,o))</'(n, m + 1, U(o4), M(i,i) ; a, /?) = 0(n + 1, m + 1, M(i,i) , M(2,i) ; a, /3) . 

Obviously, if we take (/)(n, m, "(0,0) , "(1,0) ; ck, /?) =0 and a(n; a, /?) = constant, then the last equation 
is satisfied. Thus, so far, we have proved that, independently of the rank condition ([77|l . every equation 
Q = 0, where the function Q is affine linear and possesses the D4-symmetry, admits a three-point 
symmetry with characteristic 



i?(n, m, "(0,0), "(1,0), "(-1,0); a, /^) 



^("(0,0), "(1,0); a, ^) 1 



O '"(1,0) 

("(0,0), "(1,0); a, /3) 



"(1,0) - "(-1,0) 2 
Since the function Q is symmetric, i.e. 

Q("(o,o),"(i,o),"(o,i),(i,i) ',01,13) = e(3(u(o,o),"(o,i),"(i,o),(i,i) ;/3, a) , 
we find that the characteristic of another symmetry generator is 

■K(ri, m, U(o,o),"(o,i),"(o,-i);P, a) = o«,n,o n ("(0,0), "(0,1); P, «) 

"(0,1) ~ "(0,-1) ^ 
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Now, we focus in the generic case, i.e. laukQ = 3. In this case, the characteristic of a three-point 
generahzed symmetry generator is necessarily of the form 

f h{u/Q Q),Uii Q); a, l3) 1 \ 

i?(n,m,M(o,o),U(i,o),'«(-i,o),a,/3) = a{n;a,(3) ■ '■ -/i,„ (m(o,o),W(i,o) 

V ""(1,0) — "(-1,0) ^ ■ / 

+ ^ <l>{n, m, M(o,o) , "(1,0); a, (3) , 

where the functions a(n,a,/3), 0(ri, m, U(o.o), W(i,o); Q!i /3) satisfy the determining equation (|92p . 

Actually, the function m, U(o,o)j "(1,0)! Q^i /^) should be independent of its fourth argument, 
namely U(i,o)- Indeed, differentiating the determining equation (j92p w.r.t. W(2,o) ^^'^ using (|14b[) . 
we arrive at 

'i("(o,i), ""(1,1)) {n + 1, m, U(i,o), U(2,o); a, P) + 

G'(u(i.o),'»(o,i)) <P.u,2 + l,m + l,M(i^i),M(2,i);a,/3) ^^'^^ = 0, (93) 

C'W(2,0) 

since U(2,i) depends on W(2,o) through the equation (5 = 0. Since the value W(o,i) occurs in the last 
equation through the polynomials /i(w(04), and q), ""(ci)); we set the coefficients of the 

different powers of M(o,i) equal to zero. This leads to a linear system for the derivatives of (f> appearing 
in equation (j93p . The maximal rank condition for the matrix Q ensures that the matrix of this system 
given by 

^(W(0,1):'"(1,1)) G(U(i^o)7'«(o,l)) 

'i'(w(o,i),'«(i,i)) G'(u(i^o),'«(o,i)) 
/i" (""(0,1), "(1,1)) G"'(u(i^o),"(o,i)) 



d 

where 



n,n.n=0 



5U(0,1) 



^(0,1) 

has rank 2, which implies that the function (f> is independent of its fourth argument, i.e. 

m, U(o^o),'"(i,o);a7/3) = '/'('^, w, U(o,o) ; a, /?) • 
Thus, in this case, the determining equation becomes 

r f fQ\ r , T aw hr \2 o /G(u(i,o), «(o,i))\ 
{a{n;a,P) - a{n + l;a,fJ)) /i(m(o,o), "(1,0)) »u„ o) T7 V 

\ "1^^(0,0)7 ""(1,0) J / 

+G(u(i^o),'«(o,i))<^(n,m,U(o,o);a,/3) + /i(u(o^o), '«(o,i))0(n + 1, to, M(i^o); ^) (94) 



+/i(m(o,o), "(1,0) )0(?J, ™ + l,"(o,i) ;",/?) = Q^«(i_i)0(»^ + l,"(i,i);a,/3) . 

A final comment is that the form of the function 4> is obtained in a similar manner as the one used 
to obtain the general form of the characteristic of a Lie point symmetry generator, as presented in 
Section ID The substitution of into the determining equation ([94)1 and the usage of equation Q = 
to eliminate U(i,i) in the resulting equation yield a polynomial in W(o,o)j "(1,0) ^tnd U(o,i)- Setting the 
coefficients of the different monomials equal to zero, we come up with an overdetermined linear system of 
difference equations for the unknown function a(n;a,/9) and the functions ^^(n, to; a, /3), which occur 
in the general form of the function (j). The general solution of this system delivers the three-point 
generalized symmetries and the point symmetries, as well. 
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